Topological defects in flat nanomagnets: the magnetostatic limit 



in . 
o ■ 
o ■ 

(N : 

wo: 
s . 

<\ 
m : 



o 
o 



> 
o 

00 

o 
in 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



G.-W. Chern, H. Youk, and 0. Tchernyshyov 
Department of Physics and Astronomy, The Johns Hopkins University, 
3400 N. Charles St., Baltimore, Maryland 21218 

We discuss elementary topological defects in soft magnetic nanoparticles in the thin-film geometry. 
In the limit dominated by magnetostatic forces the low-energy defects are vortices (winding number 
n = +1), cross ties (n = —1), and edge defects with n = —1/2. We obtain topological constraints 
on the possible composition of domain walls. The simplest domain wall in this regime is composed 
of two —1/2 edge defects and a vortex, in accordance with observations and numerics. 



Nanorings made out of a soft ferromagnetic mate- 
rial generate considerable interest as prospective build- 
ing blocks for nonvolatile random-access memory Q . An 
attractive feature of the ring geometry is the existence of 
two lowest-energy states in which magnetization points 
in the azimuthal direction (clockwise or anticlockwise) 
preventing the straying of magnetic field. The switch- 
ing between the two states can be accomplished by ap- 
plying the magnetic field in the plane of the ring or by 
injecting electric current. In both cases the switching is 
accomplished by nucleating a small bubble of the oppo- 
site domain and letting it expand until it occupies the 
entire ring 0. Alternatively one can view the process 
as the creation, propagation, and mutual annihilation of 
two domain walls separating the domains with clockwise 
and counterclockwise magnetization. These considera- 
tions motivate us to study the properties of domain walls 
in nanorings. 

Domain walls in magnetic nanoparticles differ substan- 
tially from domain walls in macroscopically large mag- 
nets. The main reason for that is the more prominent role 
of the surface in smaller samples. Qualitative changes are 
expected when one (or more) of the particle dimensions 
crosses a length scale characterizing the strength of fer- 
romagnetic exchange relative to that of the stray field, 
A = yA/ (iqMq, or material anisotropy, A a = y/A/K. 
Here A is the exchange constant, M$ is the equilibrium 
magnetization, and K is the anisotropy constant of the 
material |3( • The anisotropy scale A a is particularly large 
in soft materials and can be considered infinite for sub- 
micron particles. The exchange length A is in the range 
of a few nanometers. 

Previous experimental and numerical studies of do- 
main walls in submicron rings p| demonstrate that the 
ring curvature does not have a significant impact on the 
properties of domain walls. We therefore discuss the sim- 
pler geometry of a strip. Domain walls in strips were 
studied numerically by McMichael and Donahue || who 
found (at least) two different types: "transverse walls" 
in extremely thin and narrow strips and "vortex walls" 
in thicker and wider ones. Two of us have previously 
shown that the transverse walls are composite objects 
made of two elementary topological defects located at the 
opposite edges of the strip. In the limit of an extremely 



thin and (reasonably) narrow strip the energy of a do- 
main wall comes mostly from the exchange interaction, 
so that the magnet is described by the two-dimensional 
XY model with an anisotropy at the edge Q- In this 
limit, the elementary topological defects are (a) vortices 
and antivortices in the bulk of the strip carrying winding 
numbers +1 and —1, respectively, and (b) halfvortices, 
or "boundary vortices" Q , confined to the edge and car- 
rying fractional winding numbers ±1/2. A transverse 
domain wall consists of two halfvortices with opposite 
winding numbers 

In this paper and its companion |8( we discuss the 
structure and energetics of domain walls in the limit 
dominated by the magnetostatic interaction, achieved in 
strips whose width and thickness substantially exceeds 
the exhange length A. We demonstrate that the "vortex 
walls" observed in this limit are also composite objects 
containing three elementary topological defects: a vor- 
tex (winding number +1) residing between two edge de- 
fects (winding numbers —1/2). The identification of the 
elementary topological defects and implications for the 
composition of domain walls is the subject of this pa- 
per. The companion paper 8] deals with the energetics 
of composite domain walls. 

Vortex walls are stabilized when both the width and 
thickness of a strip substantially exceed the exchange 
length A |5j . In this limit the magnetostatic energy is the 
dominant contribution to the energy of a domain wall 
0, [^} and the primary force determining the shape of 
topological defects. Because the magnetostatic energy is 
a nonlocal functional of magnetization Q, energy min- 
imization is a computationaly difficult problem. There- 
fore identification of topological defects is not as straight- 
forward as in the limit dominated by exchange 6J. Fur- 
thermore, the magnetostatic energy has a large number 
of absolute minima and one must search among these so- 
lutions for one with the lowest exchange energy, making 
this a degenerate perturbation problem. 

For simplicity we will use the geometry of a thin film 
with a constant thickness t that is small in comparison to 
the width of the strip uu. In this case the shape anisotropy 
forces the magnetization M to lie in the plane of the film 
(with the possible exception of vortex cores 52). It will 
be further assumed that the magnetization depends on 
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the coordinates in the plane of the film only, 

M = M(x,y) = (M cos 0, M sin 0,0). 



(1) 



For a given configuration of magnetization M(r) its 
magnetostatic energy (no/2) J H 2 dV can be recast as 
the Coulomb energy of magnetic charges with density 
Pm(r) = -V • M = M o (sm0d x - cos9d y 0). Being 
positive definite, the magnetostatic energy has an ab- 
solute minimum of zero, which corresponds to the com- 
plete absence of magnetic charges. Thus it makes sense 
to look for low-energy states with topological defects 
among configurations with zero charge density in the 
bulk, —V • M = 0, and on the surfaces, n • M = (here n 
is the surface normal). A method for constructing such 
solutions has been discussed by van den Berg 12]. It 
yields configurations with domains of smoothly varying 
magnetization separated by discontinuities in the form of 
Neel-type domain walls. The walls acquire a finite width 
when the exchange interaction is taken into account. 

Vortex. The simplest nontrivial example is the config- 
uration with a single vortex at the origin, exp (i0(x, y)) = 
±i(x + iy)/\x + iy\. The two signs give two different val- 
ues of chirality (direction of circulation) of the vortex. In 
both cases the topological charge, or the winding num- 
ber ^3 j is +1- This solution has zero density of mag- 
netic charge and thus minimizes the magnetostatic term. 
Furthermore, it also represents a local minimum of the 
exchange energy. Taken together, these two observations 
show that the vortex is a stable configuration. Its energy 
diverges logarithmically with the system size R: 



E +1 ~ 2wAtlog (R/X) + E c 



(2) 



Antivortex. In contrast, the antivortex configura- 
tion minimizing exchange energy, exp (i0(x, y)) = ±i(x — 
iy)/\x — iy\, has a nonzero density of magnetic charge and 
thus represents a poor starting point for constructing a 
bulk topological defect with the winding number —1 in 
this limit. Minimization of the magnetostatic energy in 
this topological sector is achieved in the configuration 
known as the cross tie 14], an intersection of two 90° 
Neel walls normal to each other (Fig. QJ. The energy of 
an antivortex grows linearly with the length of the Neel 
walls L emanating from it: 



E-x ~ atL + E c 



(3) 



(the core energy is generally different from that of a vor- 
tex). The surface tension of the wall a is determined by 
the competition of exchange and magnetostatic forces. 
When the film thickness t exceeds the Neel-wall width 
(of order A), the calculation simplifies: the magnetiza- 
tion depends only on the coordinate transverse to the 
wall. The surface tension is then found by minimizing 
the total energy per unit area Q 



cr = J dx [A(d6/dx) 2 + [i Mq (cos 6 
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FIG. 1: Top to bottom: a vortex, an antivortex, and a —1/2 
edge defect in the magnetostatic limit. 



cos0 o ) 2 /2], (4) 
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FIG. 2: Determination of the topological charges at the edges 
and in the bulk. 



subject to the boundary conditions 6 (±00) = ±9q, where 
29q is the angle of spin rotation across the wall. Mini- 
mization of the total energy yields a domain wall of a 
characteristic width \\f2f sin# with surface tension 



a = sin#o — $0 cos#o) A/\. 



(5) 



In thinner films (t ^ A) the magnetostatic term is nonlo- 
cal and the Neel walls acquire long tails Q . 

Edge defects. An edge defect with the winding 
number —1/2 can be constructed by placing the core of a 
cross tie at the edge of the film, so that the magnetization 
along the edge is parallel to the boundary (Fig. HJ. As 
the core is circumvented counterclockwise the magnetiza- 
tion rotates clockwise through 7T, in agreement with the 
definition [6J. The energy of such a defect is also given 
by Eq. 

We have not been able to find any configuration that 
would contain a +1/2 edge defect and be free from mag- 
netic charges. It looks likely that the +1/2 defects carry 
a finite amount of magnetic charge and thus have a sub- 
stantially higher magnetostatic energy than the other 
three types of defects described above. This may indicate 
that, in the limit where the magnetostatic energy dom- 
inates, a +1/2 defect will decay into a vortex [n = +1) 
and an edge defect (n = —1/2). The +1/2 defects are 
stable in the exchange limit [6J, [3 ■ 

Defects and composite domain walls. The de- 
fects discussed in this paper determine the properties of 
domain walls in nanomagnetic strips. Postponing a de- 
tailed discussion to the accompanying paper [3 here we 
make two general observations that place important con- 
straints on the possible composition of a domain wall. 

First, a domain wall in a strip must contain (at least) 
one edge defect at each edge. This follows from the def- 
inition of their winding numbers [jj. Moving along the 
upper/lower edge (Fig. |5J one finds that the magnetiza- 
tion rotates through the angle — 27rni ) 2. In the presence 
of a domain wall, the edge winding numbers n\ and rii 
are half-integers. 

Second, the total topological charge of a domain wall, 
including the winding numbers of the edges and the bulk, 
must be zero. This can be seen by drawing a contour 
enclosing the domain wall (Fig. EJ) and noting that the 
total angle of rotation along that contour — 2^n\ — 2im2 
also equals 27rn, where n is the winding number in the 
bulk. Hence n + n\ + = 0. 



Thus domain walls with the smallest number of defects 
may contain (a) two edge defects with winding numbers 
+1/2 and —1/2 and no bulk defects; (b) two +1/2 edge 
defects and one antivortex; and (c) two —1/2 edge de- 
fects and one vortex. Case (a) corresponds to the trans- 
verse wall, which is indeed the lowest-energy domain wall 
in the exchange limit 0. In the opposite magnetostatic 
limit one must minimize the number of +1/2 edge defects 
(which have a high magnetostatic energy). Therefore it is 
reasonable to expect that the lowest-energy domain walls 
in this limit are of type (c). Both experimental observa- 
tions Q and numerical simulations [3 are consistent with 
this proposition. See the companion paper || for details. 

Much of the recent experimental effort in nanomag- 
netism has been devoted to the study of vortices [ill llil 
fl6j |. Given an equal (if not greater) importance of edge 
defects in determining the properties of domain walls, a 
careful examination of topological defects at the edge is 
highly desired. 

Acknowledgment. We thank C.-L. Chien, J.-G. Zhu, 
and F. Q. Zhu for discussions. This work was supported 
in part by the NSF Grant No. DMR05-20491. 



[1 
[2 
[3 
[4 

[5 

[6 
[7; 

[8 
[9 
[10; 

[11 

[12 
[13 

[14; 
[15; 
[16 



J.-G. Zhu, Y. Zheng, and G. A. Prinz, J. Appl. Phys. 87, 
6668 (2000). 

F. Q. Zhu, D. Fan, X. Zhu, J.-G. Zhu, R. C. Cammarata, 

and C.-L. Chien, Adv. Mater. 16, 2155 (2004). 

A. Hubert and R. Schaefer, Magnetic domains (Springer, 

1998). 

M. Klaui, Vaz, J. A. C. Bland, T. L. Monchesky, J. Un- 
guris, E. Bauer, S. Cherifi, S. Heun, A. Locatelli, L. J. 
Heyderman, et al., Phys. Rev. B 68, 134426 (2003). 
R. D. McMichael and M. J. Donahue, IEEE Trans. Magn. 
33, 4167 (1997). 

O. Tchernyshyov and G.-W. Chern, cond-mat/0506744. 
M. Kurzke, Calc. Var. PDE (to appear), preprint 
14/2004, Max Planck Institute for Mathematics in the 
Sciences. 

H. Youk, G.-W. Chern, K. Merit, B. Oppenheimer, and 
O. Tchernyshyov, these Proceedings; cond-mat/0508741. 
R. V. Kohn and V. V. Slastikov, Proc. Roy. Soc. London 
Ser. A 461, 143 (2004). 

A. DeSimone, R. V. Kohn, S. Miiller, and F. Otto, 
preprint 80/2004, Max Planck Institute for Mathemat- 
ics in the Sciences. 

A. Wachowiak, J. Wiebe, M. Bode, O. Pietzsch, M. Mor- 
genstern, and R. Wiesendanger, Science 298, 577 (2002). 
H. A. M. van den Berg, J. Appl. Phys. 60, 1104 (1986). 
P. M. Chaikin and T. C. Lubensky, Principles of con- 
densed matter physics (Cambridge University Press, 
Cambridge, 2000). 

M. Lohndorf, A. Wadas, H. A. M. van den Berg, and 
R. Wiesendanger, Appl. Phys. Lett. 68, 3635 (1996). 
T. Shinjo, T. Okuno, R. Hassdorf, K. Shigeto, and 
T. Ono, Science 289, 930 (2000). 

S.-B. Choe, Y. Acremann, A. Scholl, A. Bauer, A. Doran, 
J. Stohr, and H. A. Padmore, Science 304, 420 (2004). 



